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We study diffusive transport through ferromagnet - normal metal - ferromagnet (F-N-F) systems, 
with arbitrary but fixed magnetization directions of the ferromagnetic reservoirs and orientations 
of a magnetic field applied to the normal metal. For non-collinear configurations, the complex 
mixing conductance describes the transport of spins non-collinear to the magnetizations of the 
ferromagnetic reservoirs. When ImG fi 7^ 0, the total conductance of the system in the presence of 
a magnetic field can be asymmetric with respect to time reversal. The total conductance changes 
non-monotonically with the magnetic field strength for different magnetic configurations. This 
modulation of the conductance is due to the precession of the spin accumulation in the normal 
metal. The difference between the conductance of the parallel and antiparallel configurations can 
be either positive or negative as a function of the applied magnetic field. This effect should be best 
observable on Al single crystals attached to ferromagnetic electrodes by means of tunnel junctions 
or metallic contacts. 



I. INTRODUCTION 

In hybrid systems of ferromagnetic and normal metals, interesting phenomena can appear due the inter- 
play between charge and spin. The discovery of the giant magnetoresistance (GMR) effect in metallic magnetic 
multilayers, El has motivated a large number of studies on the transport properties of such systems. El The GMR is 
caused by spin dependent scattering in the system. Most studies concentrated on collinear configurations (parallel 
and antiparallel configurations). There are several papers which cover non-collinear magnetizations, both theoretically 
and experimental^ n 

Magnetoelectronic multiterminal devices reveal new physics,u but may also lead to novel applications, e.g. non- 
volatile electronics. Johnson and Silsbee investigated spin dependent effects in a 3-terminal device. u They found 
transistor effects that depend on the relative orientation of the magnetization of the ferramagnets.El More recently, a 
ferromagnetic single-electron transistor in a three terminal configuration has been realizedQ and studied theoretically.!!! 
In this case tha-source-drain current also depends on the relative orientation of the magnetizations. 

Brataas et alu give a unified semiclassical picture forj-electron and spin transport in such systems. Their formalism is 
inspired by the circuit theory of the Andreev reflect ion ,113 and is applicable to systems with non-collinear magnetization 
directions and an arbitrary number and variety of contacts between the ferromagnetic and the normal metals. However, 
the simple circuit theory of Ref. ^ only holds when the resistances of the contacts between the ferromagnetic and the 
normal metals are much higher than the resistance of the normal metal itself, thus fails when the size L of the system 
in the transport direction becomes too large. Moreover, when the size of the system is larger than the spin diffusion 
length (L 3> l s f), the presence of spin-diffusion in the normal metal requires a more complicated description with 
spatially dependent spin distribution functions. 

In the present paper, we present a study of the transport properties of simple F-N-F systems (see Fig. 1), taking 
into account different magnetizations of the ferromagnetic reservoirs and spin-diffusion in the normal metal. At 
low temperatures, spin-flip can be due to spin-orbit interactions and scattering by defects or impurities. Exchange 
scattering by paramagnetic impurities also flips the spin (see e.g. Appendix A in Ref. pd] ). The length of the normal 
metal L is assumed to be much larger than the mean- free path If, so electronic transport may be described by the 
diffusion equation. On the other hand, we allow the spin diffusion length l s f, which is the length scale on which an 
electron looses its spin in diffusive transport, to be much smaller, of the same order, or much larger than the size 
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of the system L. Under an applied bias, ferromagnetic reservoirs inject a spin-current, causing a non-equilibrium 
magnetization or "spin accumulation" in the normal metal. We are interested in the different mechanisms that reduce 
and also rotate this spin accumulation. For non-collinear configurations the physics of spin injection is more subtle 
than in the simple collinear case, since it requires generalized boundary conditions for transport through a single 
ferromagnetic-normal metal (F-N) contact .□ In general, such a contact is charaterized not only by the conventional 
spin dependent conductances G' ,G^, which describe the transport of spins collinear to the magnetization of the 
ferromagnetic reservoir, but also by the (complex) mixing conductance G^ (see Ref. ^|), that contains information 
about the transport of spins oriented perpendicular to the magnetization of the ferromagnetic reservoir. 

We are also interested in the effect of a magnetic field applied to the diffusive normal metal in arbitrary directions. 
In this case we assume that the magnetic field only couples to the spin degrees of freedom. Our approach is similar 
to the treatment of a processing magnetic field applied to a diffusive metal in Ref. [l^. 

In section II we introduce and solve the basic equations for the diffusive spin transport, showing the general 
expression for the non-equilibrium distribution function in the normal metal. In section III we discuss the boundary 
conditions of the problem. In section IV we obtain analytical expressions for the total conductance of the system in 
collinear configurations and in the absence of applied magnetic field. We also obtain analytical expressions for the 
total conductance in the case of non-collinear magnetization directions, zero magnetic field and no spin-flip scattering. 
In section V we calculate numerically the conductance in the general case. In section VI we summarize and discuss 
our results. 



II. DIFFUSIVE SPIN TRANSPORT 



When a bias is applied to our F-N-F device, a spin current is injected from the ferromagnetic reservoirs into the 
normal metal, causing a non-equilibrium magnetization or spin accumulation. For an arbitrary magnetic configuration 
of the system, the spin accumulated in the normal metal can be oriented in differents directions. If we take the spin 
quantization axis parallel to the magnetization of one of the ferromagnetic reservoirs, we need to take into account 
spins oriented perpendicular to this quantization axis, which can be described as a superposition of up (f) and down 
(|) spin states. We study a geometry invariant to translations in the lateral direction, so all quantities depend only 
on one spatial coordinate (x). The spin-polarized electron distribution is characterized by a 2 x 2 matrix in spin space 
of the form: 

/ aw m*) \ 
/» = • (i) 

V /ff(») AW / 

When the size of the system L is larger than the spin diffusion length l s f, f N (x) depends on the position. Here we 
are interested in transport under the condition If <C l s f, where If = Vp (1/r + l/r s f) 1 is the mean free path, vp is 
the Fermi velocity, r the spin-conserving scattering time and r s f the spin-flip scattering time. Both r and r s f are 
considered isotropic in momentum space. The spin diffusion length l s f is defined as l s f = *J Dr s f , where D = vplf/d, 
is the spin- independent diffusion coefficient of the normal metal (d — 1,2,3 is the dimension of the normal metal). 
So under the condition If <^ l s f, we obtain for diffusive spin transport in the steady state the following 2x2 matrix 
equations for f N (x) 



f(x) = -D 9 -^. (3) 

where 1 is the unit matrix and where the electron charge e is assumed to be equal to one. Eq. (|^) describes the 
relaxation of the spin accumulation due to spin-flip scattering, and (|^) relates the current density matrix j N (x) and 
f (x). In the case of collinear transport, our matrix equations simply reduce to 

d*f»(x) _l f s N (x)-f N s (x) 

8x 2 2 P sf { ' 

j^ x ) = -D d J^± (5) 
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where s — ("f , J,). Eqs. (ji|) and (||) have been extensively used for collinear transport in F-N multilayers in which the 
current flows perpendicular to the planes of the interfaces (CPP geometry ).liIE3 

We are also interested in the effect of an external magnetic field applied to the normal metal in an arbitrary direction. 
We know that the magnetic Zeeman energy associated with the coupling between the magnetic field and the spin of 
the electrons is given by g[iB^'B /2, where /xg is the Bohr magneton, g is the gyromagnetic ratio, a = (a x , a y , a z ) is 
the vector of Pauli matrices and B is the external magnetic field. Semiclassically, we can write for the spin dynamics 
(see e.g. Ref. |3) 



df N {x) _ i 



Then, in the steady state: 

d 2 f N (x) 



D- 



01 



9fJ>B 



{a-B) , f N {x) 



dx 2 



= - />)- 



T* (/"(»)) 



9VB 



a-B), f N (x) 



(6) 



(7) 



Using the properties of the Pauli matrices we can express the non-equilibrium distribution matrix f N (x) as: 

f N (x) = f (x)l + a-f(x) 



(8) 



where fo(x) is a scalar and f(x) — (f x (x), f y (x), f z {x)) is a three component vector. fa(x) is the particle or spin- 
independent distribution function. On the other hand, f z (x) describes the "spin polarization" on the system, and 
f x {x) and f y {x) contain information about the spins oriented perpendicular to the quantization axis. We call the three 
component vector f(x) the spin- dependent distribution function. Using (||), we separate (0) into two contributions, 
one for the spin-independent part and another for the spin- dependent part: 



d 2 fo(x) 

dx 2 
d 2 f(x) 
dx 2 







;2 



/(*) 



g^B B 
h D 



x f(x) 



(9a) 
(9b) 



The spin- independent part (Eq.(pa|)), is the conventional result for diffusive particle transport. Similar to Eq. (|^), 
the particle current density Jq (x) reads 



d 



D- 



Tr 



(/*(*)) 



dx 



-2D 



dfo(x) 
dx 



(10) 



Eqs. 



i|) and (10) express the particle current conservation 

dj^(x) . 



The general solution of Eq. 



is 



dx 



fo(x)=V + Ox. 



(11) 



Eq. (9b) describes how the spin accumulation relaxes by spin- flip scattering and by the spin precession around the 



magnetic field. This equation can be written in a general matrix form as: 

d 2 f(x) 



dx 2 

The eigengenvalues associated with the matrix A are: 



A/(:r) 



where we have introduced the vector h = gfisB/hD, which describes the "effectiveness" of the magnetic field in a 
diffusive metal. The eigenvector associated with A G is 




On the other hand, A+ and A_ have associated two complex conjugated eigenvectors v+ = V\ + iv-i and V- — v\ — iv2, 
where 



vi = 




and 



f>2 



K + ^ 




The general solution of Eq. (Bq), can then be written in terms of the eigenvalues A G , A+, A_ and vectors v , v%, V2 as: 



m = 



A v cosh(x/l s f) + B v sinh(cc// s /) 
+C [vi cosh(A) cos(F) - v 2 sinh(A) sin(F)] 
-V [vi sinh(A) sin(F) + v 2 cosh(A:) cos(r)] 
+£ [ui sinh(X) cos(F) - v 2 cosh(X) sin(F)] 
-T [vi cosh(X) sin(Y) + v 2 sinh(A") cos(F)] . 



(12) 



where 



X 



Y 



1 + Vl + a 2 x 



-1 + vt~t 



sf 



and where the dimensionless constant a = guj^Tgf = 



mechanisms. u>l = [J>b 



D 



l s f is the ratio between spin-flip and precession relaxation 

/h is the (Larmor) frequency for the spin precession around the magnetic field.0 The 

solution associated with A Q describes the relaxation of the spin accumulation due to spin-flip scattering, and the two 
complex conjugated solutions associated with A + and A_ describe the relaxation and precession of the spins due to 
the coupling with the magnetic field. The eight real constants (0, V 1 A, B, C, T>, £, T) must be determined by the 
boundary conditions. 



III. BOUNDARY CONDITIONS 



We consider two ferromagnetic reservoirs attached to a diffusive normal metal through some arbitrary contacts, as 
shown in Fig. 1. The ferromagnetic reservoirs are supposed to be large and in local equilibrium at chemical potentials 
M£,rc (A "R. denotes left and right reservoir respectively), and with energy-dependent diagonal distribution matrices in 
spin space f£ TC (e). The components of -^(e) are given by the Fermi-Dirac distribution function f FD (e, nc,iz)> and 
the direction of the magnetization in each ferromagnetic reservoir is denoted by the unit vector mc,TZ- The current 
through the system and the non-equilibrium distribution function in the normal metal are completely determined by 
the relative orientation of the magnetization directions in the ferromagnetic reservoirs, the contact conductances, the 
normal metal conductance, the spin diffusion length and the magnetic field. 

The current through an F^N contact is given in Ref. ^| in terms of the microscopic scattering matrices of the 
Landauer-Biittiker formalism.E3 According to Eq. (3) of Ref. ^, the particle current through a single contact directed 
into the normal metal can be written as 
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i c (x) = G T ti T (f F - f N (x)) u T + G L u l (f F - f N (x)) 
-G n u^f N (x)u l - (G n )*u l f N (x)u^ 



(13) 



where f N (x) and /J ^ are isotropic distribution functions, G^ and G^ are the conventional spin-dependent conduc- 
tances, which describe the transport of spins oriented in the direction of the magnetization of the adjacent ferro- 
magnetic reservoir, and G'* — ReG^ + ilmG^ is the mixing conductance, which contains information about the 
transport of spins oriented in perpendicular direction to the magnetization of the ferromagnetic reservoir. The ma- 
trices u = (l+ a ■ m) /2, and ix 1 = (l— a ■ m) /2 define the basis in which the spin-quantization axis is parallel to 
the magnetization of the ferromagnet (for details see Ref . |l7|) . Eq. ( p"3| ) relates the spin current through the contact 
i c (x) and the non-equilibrium distribution matrix f N (x) in the normal metal. Due to current conservation, Eq. ( |l3| ) 
is equal, at each contact, to the particle current per energy interval in the normal metal (see Fig. 2). The particle 
current per energy interval is related with the current density j N (x) as, i N {x) = S v DOS j (x), where S is the surface 
perpendicular to the transport direction and v DOS is the density of states of the normal metal. Using Eq. (0) i N (x) is 



i N (x)=-Sv DOS D?& 



(14) 



So we have 



for the left contact (x = + ) and 



i u (x = o + ) = i N {x = o+) 



(15a) 



(x = L~) = i u {x = L~) 



(15b) 



for the right contact (x = L ). By substituting ( |l3| ) and ( [Til ) into (15a) and ( 15b ), we obtain the boundary conditions 
for the left contact (x = + ): 



-S V DOS d 



df N {x 



dx x=o+ 
G T u T / Ar (0 + )« T + G l u l f N (0+)u l + 
G n rff N (0+)u l + (G n )*u l f N (0+)tf 
G T w T f F u 1 +G l u l f F u l 



(16a) 



and for the right contact (x — L ): 



S v DOS D 



df N {x) 



(16b) 



dx x =l- 
G ] u}j N {L-)u} + G l u l f N (L-)u l + 
G n u^f N (L~)u L + (G n )*u l f N (L-)u^ 
G^u^f^+G l u l f^u l . 

The set of parameters {(£, G l , Re G n , Im G Ti , ii\u^} i s in genera l dif ferent for each contact, but we have omitted 



the indices C and 1Z in (16a) and (|16bJ) for brevity. Eqs. ( J16a| ) and ( |l6q ) are two 2x2 matrix equations, that provide 
us a system of linear equations that determinate the eight unknown constants (O, V, A, B,C, T>, £, J-). 

From ( ^0| ) we can see that the total particle current i$ can be written in terms of one of these constants as: 



-2Sv DOS D 



df (x) 
dx 



-2DSv DOS 



-2V o1 G n O, 



where Gn = j;V DOS 1S the normal metal conductance and V Q i is the volume of the normal metal. By solving the 



system of equations (16a) and (16b), we can calculate this total particle current, iff is proportional to the difference 
between the distribution functions of the ferromagnets i$ oc (f F — f£) , times a quantity which does not depends on 
energy. From this quantity is possible to obtain the total conductance G T : 



*o W = G T (f F - f F ) 



(17) 
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where G T is in principle a function of the relative orientation of the magnetization directions in the ferromagnetic 
reservoirs, the contacts and normal metal conductances, the spin diffussion length and also of the magnetic fieldrJ 

G T =G T (rn c ,K,{G\G l ,RcG n ,lmG n } CK ,G N ,l sf ,B) . 

By studying G T for different values of these parameters, we obtain information about the physics of the spin accu- 
mulation in diffusive systems. 



IV. ANALYTICAL EXPRESSIONS. 



The properties of the contacts are parametrized by the spin-dependent conductances | G' , G*- , Re , Im G ^ } £ n - 
For collinear configurations of the ferromagnetic reservoirs (parallel and antiparallel), it is easy to obtain simple 
expressions for the conductance, which can be interpreted by simple equivalent circuits. When l s f L, there is no 
mixing between spin-up (|) and spin-down (J,) channels, and we obtain the conductance for the parallel configuration 

Gp = 7 G cfnG N + G l c G l n G N (1 

[gI + G^Gn + gIg^ {G l c + G l n )G N + G l c G l n 

and for the antiparallel configuration: 



GjG^Gn G^G^Gn 

Ap = ~ — rr t~i — n r~r ^ > 

T n ) Gn + G c G n I G c + G n \ Gn + G c G n 



On the other hand, when l s f <C L, spin-up (|) and spin-down (!) channels are completely mixed due to spin-flip 
scattering and the spin accumulation vanishes. In this limit we have: 



G° = h7r: + ^— r+ ^T , - t ■ (19) 




These expressions correspond to the simple equivalent circuits displayed in Fig. 3. Eq. ( |l8a| ) corresponds to a circuit 
in which the two spin channels are independent in the parallel configuration (Fig. 3a). Eq. ( |18b| ) corresponds to the 
anti-parallel configuration (Fig. 3b). Eq. ( fl9| ) is equivalent to a circuit with a complete mixing between spin up (f) 
and spin-down (!) channels (Fig. 3c), in which there is no difference between parallel and anti-parallel configurations. 

For symmetric contacts 

( G l = G 1 K = GT and G l c = G l n = G^j , we find analytical expressions for the conductance 
of the system for any value of L/l s f, in the parallel configuration: 



2G T G^ tanh( =£- ) + GG 



N 



Gp = 2G N = (20a) 

Gjv (4G W + G) + 2(GatG + GTG^)^ tanh(2^r) 



and in the antiparallel configuration: 



2G-T G 1 hL + GG N tanh( £-) 



G S AP = 2G N ^ j— (20b) 

G N (4Gat + G) tanh(^j) + 2(GatG + G^G^'-f- 

where G = G T + G^. 

In the limit l s f 3> L, these equations reduce to, for parallel configuration: 

s _ G T Gjy G l G N 
Gp ~ G1 + 2G N + Gl + 2G N > (21aj 



for anti-parallel configuration: 



G 



and in the limit (l s f <C L) for parallel and anti-parallel configurations: 



G ° = ra75- <22a) 

For non-collinear configurations there is no simple circuit analogy, but we can still find an analytical expression for 
the total conductance of the system as a function of the angle between the magnetizations of the different ferromagnets 
6, when l s j ^> L, at zero magnetic field (B — 0) and for symmetric contacts: 

G T (9) = 2G N (23) 

\G n \ 2 (4G 2 N {1 + cosd) + 2GG N ) + 2 Re G n G 2 N G (1 ~ cos 9) \ 



1 - 



\GU | 2 ((AG 2 N + GG N ) (1 + cos 9) + 2 {GG N + + 2 Re G^G N ({GG N + GWl) (1 - cos 6)) 



In the limit of = and 9 = tt, Eq. (B3J) simplifies to Eqs. (21a) and (|2 lbj) respectively. When the resistance of 



the normal metal is negligible compared to the contacts resistance (Gn — > oo), this reduces to 

G T m = ^(l-p 2 a tanH { 2 ) (24) 

2 \ tan 2 0/2 + M 2 /Re(r/) J 

where p = P/G = \G' — G^) /G is the polarization and r\ =J£G_^ /G is the (complex) relative mixing conductance 
Eq. (|24|) can also be obtained by means of the circuit theory. 



V. NUMERICAL RESULTS. 



The total conductance depends on the spin-dependent conductances of the contacts. We mostly set the polarization 
p = P/G = (G^ — G^) jG = 0.5 (for real metallic ferromagnets like Fe or Co, p is 0.4 and 0.35 respectively!^) , 
which corresponds to a_ ratio G^/G-^ = 3. On the other hand, the real part of the mixing conductance obeys 
ReG^ > (Cfi + G^) /2.EI The conductances of the contacts and the diffusive normal metal are considered to be of the 
same order G N - (G\G l , G n ) . 



A. Collinear and non-collinear configurations 

The total conductance depends on the magnetic configuration. We plot in Fig. 4a and Fig. 4b, G T /Gn as function of 
the relative angle between magnetizations 9, for symmetric contacts, zero magnetic field (B = 0) and in the absence 
of spin-relaxation in the normal metal (l s f L), as given by Eq. ( p3| ) for different values of ImG^ and ReG^ 
respectively. For 9 = 0°, 9 = 360° and 9 = 180° the total conductance does not depend on the mixing conductance 



and the values of G t /Gn at 9 — ,360 and 9 = 180 are given by Eq. ( |21a| ) and Eq. (21b) respectively. On 
the other hand, for non-collinear configurations, the total conductance increases with increasing mixing conductance 
(the dip become more sharp). This enhancement is due to the contributions of non-collinear spins to the transport, 
in which electrons with spins oriented in different directions than the magnetization of the adjacent ferromagnet are 
transmitted or reflected at the contact. These processes are described by the real and the imaginary part of the 
mixing conductance. 



B. Spin- flip scattering. Spin relaxation 

When spin-flip scattering is caused by spin-orbit interaction in the normal metal, the spin diffusion length l s f can 
be estimated to be equal to lf/(aZ) 2 , where a is the relativistic fine structure constant, Z is the atomic number and 
If is the mean free path (see e.g. Ref. ^). In Co/Cu multilayers, the spin diffusion length l s f is of the order of a few 
hundred angstrom (see Appendix A in Ref. |ll|) . For Al, l s f can be estimated to be of the order of a few micrometers 
for polycrystalline Al (see Ref. |l9|), or even between 10 — 70[im for Al-single crystals. In the case of very pure Na, 
r s f ~ l/isril In this case, l s f limited by spin-orbit interactions can be estimated to be of the order of 0.4 cm. In Fig. 
5a we plot, at zero magnetic fie ld a nd for symmetric contacts, the conductance of the system G T , normalized to the 



conductance G° given by Eq. (22a), as a function of L/l s f. The length of the normal metal section L is set to be 
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constant, in order to keep a constant value of Gn- When l s f 3> L the conductance of the system depends on the 
magnetic configuration. By decreasing l s f, all configurations converge to the same value of conductance G T /G° = 1. 
All configurations reach the same value of the conductance long before G T /G° — 1, since for l s f < L both contacts 
become independent and as the relative magnetic configuration is irrelevant. In Fig. 5b we plot G T /G° in the 
case of antiparallel configuration for differents values of the relative polarization P/Gn = (G^ — G^) /Gn and for 
G/Gn = (G^ + G^) /Gn constant. When l s f ^> L the configuration with large relative polarization P/Gn gives a 
small conductance and vice versa. The spin accumulation increases with increasing polarization of the ferromagnet 
and causes a reduction of the total conductance of the system. For l s f « L we also see that in each case the 
conductance approaches G° asymptotically in different ways, depending on the magnitude of the spin accumulation. 



C. Effect of the magnetic field. Precession and relaxation 

In a diffusive system the presence of an external magnetic field relaxes the spin accumulation, in addition to the 
usual precession of the spin. Semiclassically, the spin accumulation at a certain position x is the average contribution 
of the spin of all electrons. In a diffusive metal each electron diffuses along a random trajectory, while its spin precesses 
with frequency ojl around the magnetic field. Since each trajectory has a different length, the spins of the electrons at 
a certain point x are oriented in different directions, which in average relaxes the local spin polarization. The length 
scale of both relaxation and precession processes is the precession length I b — / gfisB, where D is the diffusion 

coefficient, B is the magnetic field, /is is the Bohr magneton and g is the spin gyromagnetic ratio. 

The external magnetic field may also influence the transport processes described by the mixing conductance G^ 
at the contacts. Let us consider for simplicity that Gn (G' , G^ , G'^J, and l s f — > oo. In this limit the distribution 
function of the normal metal does not depend on position. From current conservation we have: 

*£ C + ^ = (25) 



7° + T C 



(Sf-S x /) V ol (26) 



where Eq. (|25| ) corresponds to the particle current, Eq. (|26| ) corresponds to the spin current and V Q i is the volume 
of the normal metal. The current is defined to be positive when injected into the normal metal by the ferromagnetic 
reservoirs. We can re-write (^) as 



'C 



+ in = {g^L x fj V i 



where Col is the Larmor frequency vector. From this expression follows that the time scale relevant for lul is the 
escape time r esc = e 2 v DOS . V i / G contact , where G contact is the average contact conductance. r esc is the time in which 
an electron escapes from the normal metal into the ferromagnetic reservoirs. It is also the time scale relevant for the 
precession of the electrons around the magnetic field. On the other hand, if Gn ~ (G^ , G^ , G^) , r esc is of the order 
of the Thouless time td, which is the average time in which an electron passes through the diffusive normal metal. 
When Gn ~ (G' , G^ , G^) diffusion, precession and transmission or reflection at the contacts, happen on the same 
time scale. From these estimates we see that the ballistic or diffusive nature of the normal metal is not going to change 
the effect of the magnetic field on the physics at the contacts. The results obtained for Gat 3> (G^ , G^ , G^), should 
therefore be valid when G^ - (G T , G l , G n ). We now make a perturbation expansion in small magnetic fields (see 
Appendix A). To first order, the current depends on the expansion parameter B as: 

^ = S *C- 1 M Cr 1 b. (27) 
dB v ; 

where s and b are vectors associated with the spin current injected into the normal metal (see Appendix A) and 
where the matrix C describes the contacts and M the magnetic field contribution. As detailed in Appendix A, C 
has a symmetric part Sg, which only includes three of the four contact conductances, i.e., only the conductances 
G^G^ReG^ of each contact (C and 1Z) respectively. On the other hand, C has also an antisymmetric part which 

T I * 

only depends on the imaginary part of the mixing conductance of each contact ImGj- R . The matrix M, which 
describes the precession of spins due to the magnetic field, is also antisymmetric. Using the symmetry properties of 
the matrices C and M we can determine from Eq. (|27]) the symmetry properties of the total conductance of the 
system G T = io/(fr ~fn) with respect to the magnetic field B. 

When C is a symmetric matrix, dG T /dB = for small values of magnetic field. The conductance of the system is 
then symmetric with respect to a change of sign of the magnetic field {B — » —B), i.e., with respect to time-reversal. 
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On the other hand, if C is antisymmetric, dG T /dB ^ and we can expect asymmetric behavior of the conductance 
with respect to change of sign of magnetic field. 



1. Modulation of the conductance by the magnetic field. Symmetry with respect to time reversal 

In the following, we discuss the dependence of the conductance on the magnetic field. We obtain G T /G° as a 
function of L/Ib ~ L\f~B for different magnetic configurations, l s f — > oo and L constant. In Fig. 6a we plot 
G T /G° in the case of symmetric contacts, where the magnetic field is perpendicular to both magnetization directions 
B ■ rhc,Ti — 0. In this case all injected spins precess around the magnetic field. When Is 3> L, the spins injected 
from one ferromagnet are not strongly affected by the magnetic field, so they travel through the normal metal and 
reach the other ferromagnet without relaxation. As a result, the total conductance depends on the relative magnetic 
configuration. By decreasing Ib, the spin accumulation precesses and relaxes on the scale of Is- Due to the precession 
of spins, the conductance displays in general a non-monotonic behavior with L/Ib- This modulation of the conductance 
can be understood in terms of the "matching" of the spins at the contacts after precession. According to the values 
of the contact conductances for the different magnetic configurations the spins are reflected or transmitted at the 
contacts depending on its orientation. Concerning the relaxation, the configuration with more spin accumulation (in 
this case, the antiparallel configuration) is the most sensitive to the magnetic field (increases faster than the other 
ones), since there are more spins to be rotated by the magnetic field in this configuration than in others. In Fig. 
6a, the conductance of 9 = 180° (antiparallel) and 9 = 90° configurations cross the conductance for 9 — 0° (parallel) 
around L/Ib = 1 • That means that at this point the spins accumulated in these two configurations have been reduced 
to the value of spin accumulation of 9 = 0° configuration. After the point L/Ib = 1, the parallel configuration (9 = 0°) 
gives a smaller conductance and the antiparallel configuration (9 = 180 ) gives the highest conductance. As a result, 
for L/Ib > 1, the parallel configuration is more sensitive to the magnetic field than the antiparallel configuration 
(now the one which increases faster) . The relaxation of spins via the precession around the magnetic field depends 
on the amount of spin accumulation in the system. This non-mononotic behavior of the conductance is specially 
relevant between parallel and antiparallel configurations, because the difference between the conductance of both 
configurations Gp — G T AP can be modulated from positive to negative values by the external magnetic field. 



When B ■ rhc,n = 0, according to Eq.(A7) in Appendix A, B has only one component B 3 uj ~ B 3 (jfic X m-R.) (see 



Rcf. |22| ). Moreover, the spins are injected with directions along rhc and rh-ji, so the precession due to the mag netic 



field only switches the spin directions between rhc an( i As a result, the distribution function given by (A3), has 
only two components / = fiu+ $2V. In this particular case, C reduces to S^,, which is a symmetric matrix. The 
same holds for the matrix M, which reduces to its 3 x 3 upper box, which only includes B3 (see Appendix A). As C 
reduces to S^,, we expect dG T /dB — 0. Fig. 6b shows the dependence of G T /G° on B/Bry for different magnetic 
configurations, where Bp = 2H/gfiBT D is the scale of magnetic fields relevant for precession in a diffusive medium. 
As expected, all configurations are symmetric with respect to a change of sign in magnetic field (B — ► — B). 



2. Modulation of the conductance by the magnetic field. Asymmetric properties with respect to time reversal 

Now we want to investigate the role of ImG^. To this end, the magnetic field is assumed to be oriented perpen- 
dicular to both magnetizations when the system is in collinear configurations, and paralle l to the direction of one 



of the magnetizations when the system is in the 9 = 90° configuration. According to Eq.(A7) in Appendix A, for 
9 = 90°, the magnetic field is along B\u + B2V, and as a result from the injection and precession, there are spins 
in the three directions / = f\u + /2V + /3W, i.e., the precession of spins around the magnetic field induces spins 
along the perpendicular direction (rhc x wr) to the injection orientations rhc and rhfi. C is then an antisymmetric 
matrix, due to the contributions of the terms which includes ImG^. So for ImG^ 7^ 0, dG T /dB 7^ 0, which means 
asymmetric behavior of the conductance with respect to time reversal. On the other hand, if we put ImG^ = 0, C 
is symmetric and dG T /dB = 0. 

In Fig. 7a we obtain G T /G° vs L/Ib, for the same set of parameters as in Fig. 6a. For parallel and antiparallel 
configurations, the results are not modified compared to Fig. 6a. However, for 9 = 90° the relative conductance 
G T /G° does not approach unity asymptotically. In this configuration there are some injected spins, which are parallel 
to the magnetic field and which do not precess at all. So this part of the spin accumulation remains in the system 
and does not relax irrespective of the values of the magnetic field. More interesting is the appearence of a dip in the 
conductance for small values of the magnetic field. If we repeat the calculation of G T /G° vs L/Ib for the same set 
of parameters except for ImG^/Giv = 0, we see that the dip disappears (Fig. 7b), so according to our discussion, 
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it is related with asymmetric properties of the conductance. In Fig. 7c, we plot G T /G° vs B/Bo- As we expect, 
9 = 90° configuration presents asymmetric behavior respect time reversal, whereas both parallel and antiparallel 
configurations remain symmetric (Fig. 7c). In particular the 9 = 90° conductance is antisymmetric with respect to 
time reversal, for small values of magnetic field. 

From this discussion, we understand that the real part of the mixing conductance describes processes at the 
contacts in which spins perpendicular to the magnetization direction, are transmitted or reflected obeying time- 
reversal symmetry. On the other hand, the imaginary part of the mixing conductance describes processes in which 
the spins precess around the magnetization vector of the ferromagnet. As a result of the precession, the orientation 
of the spin changes. The latter processes are antisymmetric with respect to time-reversal. 

3. Supression of the magnetic field effects by spin-flip scattering 

Spin-flip scattering causes relaxation of the spin accumulation in the normal metal and as a result, supression of 
the spin-dependent properties on the system. Now we want to investigate how the spin-flip affect the magnetic field 
effects show above. The existence of spin-flip scattering reduces l s f. If l s f ^> L, there is no strong spin-flip scattering 
in the system and it is possible to observe spin-dependent effects. On the other hand, if l s f <C L, the injected spins 
relax very fast due to spin-flip processes and no spin-dependent effects can be observed. In particular in Fig. 8 we 
show how the dip of 9 — 90° configuration from Figs. 7a and 7c, is suppressed by spin-flip scattering on the system. 
Also by decreasing l s f, G T /G° increases for constant magnetic field, to the value 1. That simply means, that the spin 
accumulation relaxes due to spin-flip scattering, as is expected. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper, the normal metal in our F-N-F device is considered three dimensional, but can also be two dimensional 
(2D), e.g., a two dimensional electron gas (2DEG) attached to ferromagnetic reservoirs, or even one dimensional 
(ID), if the normal metal isr-sti quantum wire or a carbon nanotube.Ej In this case electron-electron interaction 
should be taken into account.cil The non-magnetic material can also be a semiconductor, as shown in recent spin- 
injection experimentsEj In the case of a 2DEG attached to metallic ferromagnets, the large difference between the 
conductivities of the 2DEG and the ferromagnetic reservoirs suppresses the spin-injection via metallic contacts. For a 
significant spia-injection into the 2DEG, tunnel contacts, a semiconductor ferromagnet or a half-metallic ferromagnet 
are required E3 

In this paper we have shown how the spin-dependent transport through a F-N-F double heterojunction can be 
described in terms of the spin dependent conductances of the contacts (G^ , G* , G^) , the magnetization direction 
fh of the ferromagnetic reservoirs, and the normal metal conductance Gj\r- The dependence of the conductance on 
the relative angle between the magnetizations of the different ferromagnets is affected by the mixing conductance 
G*^. For non-collinear transport between the ferromagnetic reservoirs, G^ = Re G^ + i Im G^ describes transport of 
spins perpendicular to the magnetization direction of the ferromagnets. These processes enhaces the conductance fc* 
non-collinear configurations, which may be used in multiterminal devices for modulation of the transport properties.^ 
This modulation could be useful for future applications as spin-dependent transistors. We find that spin injection 
can be symmetric and antisymmetric with respect to time-reversal. The symmetric processes are described by Re G^ 
and the antisymmetric ones are described by ImG^. It is interesting to observe that the antisymmetric processes 
described by ImG^ correspond to spin precession around the magnetization vector of the ferromagnet which couples 
to an external magnetic field. 

In a diffusive system, an applied magnetic field produces both precession and relaxation of the spin accumulation. 
The conductance displays a non-monotonic behavior on the scale of the precession length Is, which is the distance 
for the precession of the spin around the magnetic field in the normal metal. Due to this modulation, the difference 
between the conductances of the parallel and antiparallel configurations Gp — G T AP can be positive and negative as a 
function of the magnetic field. A possible candidate to observe this effect is Al, which has a large l s f and which can be 
coupled to ferromagnetic reservoirs (e.g., Fe, CoFe, NiFe, Co,..) via metallic junctions or also AI2O3 tunnel junctions. 
Let us estimate the values of magnetic fields for Al-single crystal associated with the points L/Ib — 0.5 and L/Ib = 2 
of Fig. 6a, where Gp — G T AP is positive and negative respectively.If the length of the system is L = 10^m, which is 
comparable with the spin diffusion length (L ~ l s f = 10 — 70^m)Jj we obtain for L/Ib = 0.5 : B^ ~ 0.01T and for 
L/l B = 2 : BH - 0.1T. Ifi = lfj.m(L < l sf ), we obtain for L/l B = 0.5 : - 0.1T and for L/l B = 2 : B^ ~ IT. 
In both cases we see that the values of magnetic field in the case of Al-single crystal are quite reasonable and also 
that the change from positive to negative Gp — G T AP can be achieved by an increase of the magnetic field by one order 
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of magnitude. The same estimate for polycrystalline Al gives us for L *~ l s f = 0.7/xm, the values of flW - IT 
and ~ 16T respectively. These fields are much higher than the typical switching field for a ferromagnet, so 

polycrystalline Al does not appear to be a good candidate. For very pure Na, if L ~ l s f — 0.4cm, the corresponding 
values of magnetic fields are BW ~ //T and fit"! - 10/xT respectively. This modulation of G T P - G T AP by a magnetic 
field can also be explored in semiconductors (SC) 2DEG, as e.g. GaAs and InAs. When Is ~ l s f the following 
expression holds for the magnetic field corresponding to L/Ib = 1 : B = (t s / g) -1 = 2.27 • 10~ n (r s / ff) - , which 
depends on the spin-flip time r s / and on the gyromagnetic ratio g of the semiconductor material. For SC, g depends 
strongly on the material (e.g., g GaAs = — 0A,g InAs = 15.0) , sodepending on the values of r s /, one can obtain the 
corresponding values of magnetic field. Kikkawa and AwschalomEa report r s f ~ 10~ 7 s in n-type GaAs system, but this 
values corresponds to spin lifetimes of optically pumped carriers, and not to the usual carriers relevant for transport. 
The corresponding value for the magnetic field for this case is B GaAs ~ 5 10~ 4 T. On the other hand, we are not 
aware of reliable values of r s j for transport in these systems. In conclusion, from our estimates of the relevant values 
of magnetic fields, Al-single crystals with ferromagnetic contacts are good candidates to test our predictions and 
possibly lead to the discovery of other new physical phenomena of spin-transport. 
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APPENDIX A: PERTURBATION EXPANSION IN SMALL MAGNETIC FIELDS 



Eqs. (|25| ) and (|26|) can be written as follows: 



+ + % if' ™c) + % (f ■ ran) = %fl + %/| (Al) 



2 2 ' 2 V J 2 



°j- - Re C#J (/• rh c ) mc + Rc G% f + Im G% [f x m c 

_ ReG^ (/ • m n ) m n + Re G# f + Im g£ (/ x m K ) (A2) 

Now we expand / into a convenient basis of the vectors mc, m-ji, and mc x mn as 

f = fiu + f 2 v + hQ (A3) 

where 

- "^+^ ( A4) 

(A5) 
(A6) 



^2(1 


+ m) 


m c - 


■mn 


\/2(l 


+ m) 


mc x 


m n 


Vi- 


m? 



and where m = rhc ' m iz = cos 9. We can also express the magnetic field in this basis as 

B = Biu + B 2 v + B 3 w. (A7) 



In terms of this expansion, we can combine (Al) and ( |A2| ) into a compact matrix form as 

f C + M) a = b (A8a) 
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where 



, ■ 0\ 

^(imG^-ImGy) 

(im G 1 C L - Im G^ ) (lm G [ L + Im G# ) Re G j. 1 + Re G. 



Sc 

(3 x 3) 



-, 



-.Tl 
T n 



(A8b) 



S 6 = 



l+m Pc+Ptz 
2 2 



m Pc-Ptz 



^ fl+^ Pc+Ptz (Gc+G^jl+m) | (Re G^+R.c G^)(l-m) ^ G r. — Gt?. Rc , R c c] + ) ( ^V"^ ) 

^ (G^-R C G^+RcG^)(^i^) (G.+G^d-^) + (ReG^ + RcGn)(l +m) ^ 



and where 



M = V ol 
































Si 
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-Si 






(A8c) 



(A8d) 



a = 



( fo 
fi 
h 



(A8e) 



/ _Gcfc+ G n fi \ 



V 



^(P£ /£ - Pn fi) 







(A8f) 



By a perturbation expansion in small magnetic fields, we may study how the magnetic field is coupled with the physics 
at the contacts. To zeroth order in magnetic field we simply have 



a(°) - C- 1 b. 

To first order 

a« = (C- X M C- 1 ) b. 
The total particle current in the system is given by 

*o = i°c - in = s-a + Gc fl + G n fi 

where 

/ _G C + G n \ 

^(-Pc + Pn) 



(A9) 



(A10) 



(All) 



^(Pc + Pn) 



(A12) 
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The dependence of the current on the expansion parameter B is given by 

d in _ _ 
-dB =S - a 

which to first order reduces to 

^- = s-a^ =sC- 1 MC- 1 b. (A13) 
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Figure Captions 

Fig. 1: Ferromagnetic reservoirs attached to diffusive normal metal through arbitrary contacts. Arbitrary but fixed 
magnetization directions of the ferromagnetic reservoirs and orientations of a magnetic field applied to the normal 
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metal are taken into account. The contacts are described by the spin-dependent conductances G^,G^,G^ and the 
normal metal is characterized by the normal metal conductance Gn- The ferromagnetic reservoirs are supposed to 
be large enough and in local equilibrium. 

Fig. 2: Current conservation imposes boundary conditions on the system. The current through the contacts i c (x) 
is equal to the current into the normal metal i N (x) at each contact. i c (x) depends on the contact conductances and 
on the direction of the magnetization of the adjacent ferromagnet reservoir and i N (x) is the current for the normal 
metal. 

Fig. 3: Equivalent circuits for parallel and antiparallel configurations in the limits l s f 3> L and l s f <C L. (a) 
and (b) correspond to the parallel and anti-parallel configuration respectively, when l s f 3> L. In this limit, the two 
spin channels are independent and there is no mixing between them, (c) parallel and anti-parallel configurations 
when l s f <C L. In this case, there is complete mixing between spin-up (|) and spin-down (J.) channels and the spin 
accumulation vanishes. 

Fig. 4: Mixing conductance: Dependence of G T /Gn on the relative angle 9 between the magnetizations of 
the ferromagnetic reservoirs, for symmetric contacts, zero magnetic field and in the absence of spin-flip scattering, 
(a) The following set of parameters is chosen: G^ /Gn = 1.0, G l /G N = 0.3, Re G n /G N = 0.7 and lmG n /G N 
takes values 0.0, 0.5, 1.0, 2.0 and 10.0 corresponding to the different plotted lines, (b) In this case, G^ /Gn = 
1.0, G l /G N = 0.3,ImG Ti /GAr = 0.0 and ReG tl /G N changes with values 0.7, 1.0, 2.0 and 10.0. According to the 
condition RcG n > (G T + G l ) /2, Re G n /G N cannot be, smaller than 0.65. 

Fig. 5: Effect of spin-flip scattering on the system: For symmetric contacts and zero magnetic field, (a) 
G T normalized to G° as a function of L/l s f, for the following set of parameters: G^ /Gn = 1-0, G^ /Gn = 0.3, 
ReG^/G N - 0.7, lmG^/G N = 0.0. (b) G f /G° versus L/l s f in the case of antiparallel configuration, for different 
values of the relative polarization P/Gn = 0.1, 0.7, 1.1 and for G/Gn = 1-3 constant. 

Fig. 6: Magnetic field dependence in the absence of spin-flip scattering: We consider symmetric contacts 
and the following set of parameters: G^/Gn = 1.0, G L /G N = 0.3, ReG n /G N = 0.7, lmG n /G N = 0.5. Moreover 
the magnetic field is always perpendicular to both magnetizations directions B ■ rn^n — 0- (a) G T /G° as a function 
of L/l B . (b) G T /G° versus B/B D . 

Fig. 7: Magnetic field dependence in the absence of spin-flip scattering: We consider B ■ rhc = 0, 
B ■ m£' 180 ° = and B \\ (or B ■ m n = 0, B ■ m°' lm ° = and B\\ mf°). (a) G T /G° vs L/l B , for symmetric 

contacts and the following set of parameters: G t /Gat = 1.0, G l /G N = 0.3, ReG^/G^ = 0.7, ImG Ti /GAr = 0.5. (b) 
Same as (a) but for lmG n /G N = 0.0. (c) G T /G° versus B/B D , for G t /Gat = 1.0, G l /G N = 0.3, ReG n /G N = 0.7, 
Im G^/Gn = 0.5. 

Fig. 8: Magnetic field dependence and spin-flip scattering: Conductance for 6 = 90 configuration of Fig. 
7a, for different ratios l sf /L = 0.1, 0.3, 1, 3, 10, 100. 
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